Isospin Asymmetry in the Pseudospin Dynamical Symmetry 
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Pseudospin symmetry in nuclei is investigated considering the Dirac equation with a Lorentz 
structured Woods-Saxon potential. The isospin correlation of the energy splittings of pseudospin 
partners with the nuclear potential parameters is studied. We show that, in an isotopic chain, the 
pseudospin symmetry is better realized for neutrons than for protons. This behavior comes from 
balance effects among the central nuclear potential parameters. In general, we found an isospin 
asymmetry of the nuclear pseudospin interaction, opposed to the nuclear spin-orbit interaction 
which is quasi isospin symmetric. 
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In some heavy nuclei a quasi-degeneracy is observed 
between single-nucleon states with quantum numbers (n, 

j = i + 1/2) and {n - 1, i + 2, j = £ -fr 3/2) where 
n, £, and j are the radial, the orbital, and the total an- 
gular momentum quantum numbers, respectively. This 
doublet structure is better expressed using a "pseudo" 
orbital angular momentum quantum number, d = £ -\- 1, 
and a "pseudo" spin quantum number, s = 1/2. For 
example, for [ns]^/2i ~ l)'^3/2] one has £ — I, for 
[np3/27 {n — l)/5/2] one has £ = 2, etc. Exact pseudospin 
symmetry means degeneracy of doublets whose angular 
momentum quantum numbers are j = £ ± s. This sym- 
metry in nuclei was first reported about 30 years ago [Q , 
but only recently has its origin become a topic of intense 
theoretical research. 

In recent papers ||i §i |§7 1§1 possible underlying mech- 
anisms to generate such symmetry have been discussed. 
We briefly review the main points of these studies. 

Blokhin et al. performed a helicity unitary transfor- 
mation in a nonrclativistic single-particle Hamiltonian. 
They showed that the transformed radial wave functions 
have a different asymptotic behavior, implying that the 
helicity transformed mean field acquires a more diffuse 
surface. Application of the helicity operator to the non- 
rclativistic single-particle wave function maps the normal 
state {l,s) onto the "pseudo" state (l,s), while keeping 
all other global symmetries |^ . The same kind of unitary 
transformation was also considered earlier by Bahri et 
al. to discuss the pseudospin symmetry in the nonrcl- 
ativistic harmonic oscillator. They showed that a partic- 
ular condition between the coefficients of spin-orbit and 
orbit-orbit terms, required to have a pseudospin symme- 
try in that non-relativistic single particle Hamiltonian, 
was consistent with relativistic mean-field (RMF) esti- 
mates. 

Ginocchio for the first time, identified the pseudo- 



spin symmetry as a symmetry of the Dirac Hamiltonian. 
He pointed out that the pseudo-orbital angular momen- 
tum is just the orbital angular momentum of the lower 
component of the Dirac spinor. Thus, the pseudospin 
symmetry started to be regarded and understood in a 
relativistic way. He also showed that the pseudospin sym- 
metry would be exact if the attractive scalar, S, and the 
repulsive vector, V, components of a Lorentz structured 
potential were equal in magnitude: 5 + ^ = 0. Under 
this condition, the pseudospin symmetry was identified 
as a SU(2) symmetry of the Dirac Hamiltonian [Q. 

In RMF models, nuclear saturation is explained by a 
cancellation between a large scalar (S) and a large vec- 
tor (V) fields § . Typical values for these fields in heavy 
nuclei are of the order of a few hundred MeV (with op- 
posite signs), their sum providing a binding potential of 
about 60 MeV at the nucleus center. Therefore, a natu- 
ral claim of Ginocchio was to regard the quasi-degenerate 
pseudospin doublets in nuclei as arising from the near 
equality in magnitude of the attractive scalar, 5*, and the 
repulsive vector, V, relativistic mean-fields, S ~ —V, in 
which the nucleons move. When the doublets are degen- 
erate, the shape of the lower components of the Dirac 
spinor for the two states in the doublet is the same Q . 

More recently, Meng et al. |^ showed that pseudo- 
spin symmetry is exact when dE/dr ~ where S = 
S -\- V. They also related the onset of the pseudospin 
symmetry to a competition between the centrifugal bar- 
rier and the pseudospin-orbit potential. 

Actually, since in nuclei S{r), V{r) when r — > oo, 
the conditions E = and cE/dr = are equivalent. 
However, these conditions cannot be realized in nuclei, 
because they imply that no bound states exist Q . Never- 
theless, as we will show, there is a correlation between the 
pseudospin splitting and the depth of E, its surface dif- 
fuseness and its radius. Depending on the actual choice of 
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these parameters, fitted to describe a nucleus, the pseu- 
dospin symmetry may show up. Thus, one may argue 
that pseudospin symmetry is a dynamical symmetry in 
nuclei, in the sense of Arima's definition of a dynamical 
symmetry [^| : (i) a symmetry of the Hamiltonian which 
is not geometrical in nature; or (ii) a ordered breaking 
symmetry from dynamical reasons. This is consistent 
with the findings of Ref . ||^ , associating the pseudospin 
symmetry with a particular relation between the coeffi- 
cients of spin-orbit and orbit-orbit terms in RMF models. 

To establish that correlation, we perform a model cal- 
culation with a Lorentz structured potential of Woods- 
Saxon type in the Dirac equation. The scalar and vector 
components of this potential are the mean-field central 
nuclear potentials. The form of the potential is 

U(r) = ^ fl) 

^ ' 1 + exp[(r - i?)/a] ' ^ ' 

where U (r) stands for either the vector or the scalar po- 
tential. Although this is not a full self-consistent rela- 
tivistic potential, the use of as a nuclear mean-field 
is realistic enough to be applied to nuclei and enables us 
to study the splittings of pseudospin partners not only 
as a function of the mentioned depth, C/q, but also of the 
diffuseness, a, and the radius, R. 

In nuclear relativistic mean-field theory, each nuclcon 
is described by a Dirac Hamiltonian of a particle of mass 
m in an external scalar, S, and vector, V, potentials: 

H = a-p + p{m + S) + V , (2) 

where a and (3 are the usual Dirac matrices. The Dirac 
Hamiltonian is invariant under an SU('2) transformation 
for two cases: S = V and S ^ ~V 

We define — V — S and denote the upper and lower 
components of the Dirac spinor by = Assum- 
ing that S and V are radial functions, the Dirac equation 
can be decoupled into two equations for the lower and 
upper components, respectively: 

v^f.+ , (fI+^—^f) + 

L — m — L, \ r J 

{E + rn- A){E-m-T.)Fi=0 (3) 

E + m — A \ r J 

[E + m- A){E -m-T.)Gi = Q, (4) 

where the primes denote derivatives with respect to 
r. The spinors have been factorized in radial 
and angular parts: = iGi(r) $^(6', (^) and ~ 

—Fi{r) <3'~(0, (/)), with i standing for the quantum num- 
bers of the single particle state. The property cr • L^f = 
— (1 ± Ki)^f was used, cr being the Pauh matrices. 

Before we present our results, let us briefly discuss how 
the pseudospin symmetry gets broken. The commutator 
of the SU(2) generators of pseu dosp in symmetry with the 
Hamiltonian (|) is given by @, [uf 

[H,S.,]=^^Y'^ S) , (5) 
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FIG. 1: Pseudospin energy splittings for the neutron pseu- 
dospin partners (2/5/2,3^3/2) and (1/19/2,2/7/2) in '^"^Pb as 
a function of R for several values of a. 

where Si = ^Si^ = ^^Pi - Si and Si = crj/2. 
The breaking of the pseudospin symmetry can thus be 
related to the commutator [S,Si]. Requiring [S, Si] = 
is equivalent to the previous condition d'E/dr = 0, when 
S is a radial function. 

Now we turn to the presentation of our results obtained 
using the Woods-Saxon potentials (Q) in eqs. (||) and 
(^). There are altogether six parameters for S and A, 
namely the central depths, Sq Aq, two radii and two 
diffuseness parameters. We observed that the pseudospin 
splitting is not sensitive to R and a of the A potential, 
and, accordingly, set the same radius, i?, and surface 
diffuseness, a, for both potentials. We first fitted these 
parameters to the neutron spectra of ^°*Pb, obtaining a 
good agreement with the results of Ref. Q for the same 
set of pseudospin doublets: (liii/2, 2,99/2), (2/5/2, 8^3/2) 
and (I/19/2, 2/7/2). 

Keeping Eq and Aq fixed, we varied a and R in or- 
der to see how the energy splittings of the pseudospin 
doublets change with the surface diffuseness and the ra- 
dius. This dependence is presented in Fig. 1. As a 
increases, the splittings of the pseudospin doublets de- 
crease. On the contrary, as R increases the pseudospin 
splittings increase. The same figure shows that the split- 
tings are more dependent on R than on a. If a increases 
and/or R decreases enough, the pseudospin doublet part- 
ners cross each other, inverting the sign of the energy 
splitting, and are driven apart if a (i?) further increases 
(decreases). This inversion of pseudospin partner split- 
tings as a increases, 1+1/2 ^ 1-1/2 changing to 
^n-i 1+1/2 /-1/2' observed experimentally and 
was also found in P, p^. 

We now keep a, R and Aq fixed to study the sensi- 
tiveness of the pseudospin doublets with Eq. The results 
are presented in Fig. 2 and one observes that, as |So| 
decreases, the splitting also decreases. This is in accor- 
dance with Ginocchio predictions for pseudospin symme- 
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FIG. 2: Pseudospin energy splitting in ^"^^Pb as a function of 
the depth of the E potential for the same neutron pseudospin 
partners (similar to Fig. 1). The vertical line corresponds 
to the fitted Woods-Saxon parameters for ^°*Pb: Eo = —66 
MeV, Ao = 650 MeV, 7? = 7 fm and a = 0.6 fm. 

try breaking due to the finiteness of the S mean field. We 
also found that, for deeper levels, an inversion of pseu- 
dospin energy splitting occurs for sufficiently low |So|- 
Comparing Figs. 1 and 2, we see that the splittings of 
pseudospin partners arc more sensitive to the radius and 
diffuscness than to the central depth of the potential. 

We note that the radius dependence is relevant for 
comparing different nuclei, especially certain isotopes 
However, we found an interesting correlation be- 
tween I So I and R. Varying |So| and R but keeping a 
and the product |So|i?^ fixed, the pseudospin doublet 
splittings remain almost constant. Therefore, only two 
parameters, out of (Sq, R, a), are free. 

The variation of the other free parameter, Aq, does 
not qualitatively change the splittings. Since the A po- 
tential is related to the effective mass of the nucleons, 
the basic effect is to slightly change the nucleon separa- 
tion energies, and especially the spin-orbit coupling, as 
discussed below. We have also performed the calculation 
for calcium isotopes, and found the same dependence of 
the pseudospin splitting with i?, a and Sq. 

We apply now the systematics to some nuclei stud- 
ied in the literature. Recently, the pseudospin symmetry 
in Zr and Sn isotopes was investigated as a function of 
the number of nucleons |6|. The form of the S poten- 
tial in dependence of the radial distance for such nuclei, 
as A increases, is given in Ref. jl^, starting with ^""Sn 
and going up to ^^°Sn. In [jl^ it is shown that the cen- 
tral neutron depth, Sq, of Sn isotopes varies from about 
-65 MeV to -54 McV, R from 5.6 fm to 6.8 fm and the 
surface diffuscness increases. For the proton, the cen- 
tral depth, So, varies from -50 MeV to -60 MeV, R 
from 5 fm to 6.5 fm, while the surface diffuscness also 
increases, now slightly more than in the neutron case. 
Using (|^), good fits to the neutron central potentials of 
the Sn isotopes were obtained, with Sg from —68.6 MeV 



to —55.7 MeV, R from 5.6 fm to 6.7 fm and a from 0.6 fm 
to 0.76 fm, in agreement with Ref. 

This information on S for Sn isotopes allows us to an- 
alyze the behavior for neutron pseudospin doublets going 
from A = 100 to A = 170. As A increases, the central 
depth I So I decreases and the surface diffuseness increases, 
both effects favoring the pseudospin symmetry, as shown 
in Figs. 1 and 2. However, the radius increase with A can 
partially offset those effects. Since the values of |So|i?^ 
are roughly constant for neutrons [l^ , the correlation 
between these two parameters, mentioned above, implies 
that the effects of increasing R and decreasing |So| in 
the neutron central potential, when A increases, balance 
each other. Thus, the dominant effect comes from the 
increasing a , slightly favoring the pseudospin symmetry. 

We extended our analysis to the proton spectra. For 
this case, as mentioned above, |So| increases as A in- 
creases for Sn isotopes therefore disfavoring pseu- 
dospin symmetry. Moreover, R increases, also disfavoring 
the pseudospin splitting. Hence, at least for the proton 
spectra, and in opposition to the neutron spectra, we do 
not have any balance effect coming from So and R de- 
pendences. Thus, we have a competition between the ef- 
fects of surface diffuseness (favoring the symmetry), and 
of depth and radius of the central potential (both disfa- 
voring the symmetry). From our systematics described 
before we expect that effects coming from the depth and 
radius of the central potential override the surface diffuse- 
ness dependence. Therefore, although the proton spectra 
along the Sn isotopic chain was not presented in Ref. , 
our systematics predicts that the splitting increases for 
the proton pseudospin partners. The different behavior 
of the splitting for protons and neutrons can be inferred 
by looking at the value of |So|i?^ for the two cases as A 
changes: the change is considerably higher for protons. 

The reason for the different behavior of the central 
potential depth for neutrons and protons as a function 
of A, lies on the p meson interaction, which is repul- 
sive for neutrons and attractive for protons, as explained 
by mean-field model calculations 1 14 . The inclusion of 



this interaction, which is important in asymmetric nuclei, 
changes the vector part of the nuclear potential V to 
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with + and — signs for protons and neutrons respectively; 
Po is the time component of the p field, which is propor- 
tional to the number of protons minus the number of 
neutrons, and K; comes from the vcctor-isoscalar lo me- 
son. The increase of the parameters R and a with A can 
be traced back to the known nuclear radius A^^^ depen- 
dence and to the excess of neutrons on the surface (the 
neutron skin effect). Hence, we may conclude that, for 
a given nucleus, the parameters So, a and R for protons 
and neutrons are different. Then, an isospin asymmetry 
in the pseudospin interaction is expected to take place 
in agreement with the analysis of the mean-field nuclear 
parameters for protons and neutrons done in . In par- 
ticular, since po is negative for heavy nuclei, the vector 
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potential is bigger (and, thus, |S| smaller) for neu- 
trons than for protons. From our previous analysis, for a 
neutron rich nucleus, the pseudospin symmetry for neu- 
tron spectrum is favored, in agreement with the results 
presented by Lalazissis et al. p5[ . 

The systematics discussed here seems to be quite gen- 
eral and a comment on how it affects the spin-orbit split- 
ting is pertinent. Variations of central depth and sur- 
face diffuseness of the nuclear potential do not change as 
much spin-orbit splittings as they do for the pseudospin- 
orbit splittings. The reason is that to change substan- 
tially spin-orbit splitting one needs a significant relativis- 
tic content, i.e. a large lower component in the Dirac 
spinor The spin-orbit splitting is completely corre- 
lated with the nucleon effective mass From Eq. (jj), 
we see that it will depend strongly on A. In RMF models 
this potential carries a quite large scale when compared 
with E (Ao is around 650-750 MeV, while |Eo| is around 
50-70 MeV). Therefore, the p meson potential, V^, in the 
range 4-8 MeV, becomes irrelevant compared to A. This 
was verified by recent numerical calculations ^ , 

showing that the difference between the values of A for 
protons and neutrons (and, therefore, between the cor- 
responding values of the spin-orbit term) is very small. 
This explains why the spin-orbit interaction is roughly 
isospin symmetric. 



On the contrary, the pseudospin splitting can change 
when there is little relativistic content, i.e. a small lower 
component. This may explain why the nonrelativistic 
treatment of Ref. to analyze the origin of the pseu- 
dospin symmetry works well in explaining the small pseu- 
dospin splitting. The pseudospin-orbit interaction de- 
pends on E, as one can see from Eq. (^). Since Vp cannot 
be neglected compared to E, the different sign contribu- 
tion of Vp in for protons and neutrons has a relevant 
effect on E, leading to the isospin asymmetry in the pseu- 
dospin interaction. 

The systematics observed in our model calculation 
helps us to understand the origin of the pseudospin 
symmetry and explains the quasi-degeneracy of single- 
particle states occurring in the spectra of some finite 
nuclei. The new radioactive nuclear beam facilities will 
provide more data to which these results can be applied. 
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